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HEAD TO BASE TRANSFER FUNCTION CHARACTERISTICS OF A 
CYLINDRICAL TANK PARTLY FILLED WITH LIQUID 

by Carl F. Lorenzo, William G. Costakis, and EvaT. Jun 
Lewis Research Center 

SUMMARY 

This study presents an analysis of the transfer function of head motion to base mo- 
tion for a thin-walled cylindrical tank partially filled with a liquid. This transfer func- 
tion is part of the overall structural description for the analysis of POGO type instabil- 
ities in rocket vehicles. 

The analysis utilized the double Laplace transform to achieve a closed form approx- 
imate solution. An experimental verification of the analysis was performed using a 
cylindrical plastic tank containing water. 

The results clearly demonstrated the importance of considering fluid resonance ef- 
fects on head to base response. Indeed, for at least two f luid heights, the fluid mode ef- 
fects dominated the mechanical modes for this transfer function. 


INTRODUCTION 

The dynamic instability of launch vehicles due to coupling of the yehicle structure 
(in the longitudinal mode) and the vehicle propulsion system continues to be a problem. 
Prediction and solution of this instability known as POGO requires dynamic descriptions 
of the vehicle structure, tankage, feed system, and engine and their intercouplings. 

Significant analyses have been made into the overall stability problem (refs. 1 to 4). 
These studies generally assume that vehicle structural modes are available from a com- 
plete structural analysis or from vibration testing. A number of related studies dealing 
with the fluid-tank dynamic characteristics may be found in references 5 to 7. In par- 
ticular, references 5 and 6 consider the tank structure and fluid combination as an equi- 
valent single lump mode. Also the transfer function characteristics of the distributed 
fluid alone have been considered (ref. 7). 



The fundamental structural description required for POGO analysis is the transfer 
function of tank base pressure (feedline inlet pressure) to engine gimble plane motion 
(fig. 1). For the vehicle upper tank this requires upper tank base motion as an input, 
from which tank base pressure can be determined. This in turn requires the transfer 
function from lower tank head motion to gimble plane motion, that is, head to base re- 
sponse for a cylindrical tank. Consideration of this type response has been made for a 
cylindrical shell connecting two masses (without fluid) (ref. 8) . Also the excellent 
study of Kana and Chu (ref. 9) has numerically determined the head response with con- 
tained fluid. 




Structure dynamics 


(b) Block diagram. 


Figure 1. - Simplified block diagram of idealized vehicle. 


When the overall problem is considered in its simplest terms, the block diagram 
for POGO analysis of an idealized vehicle is shown in figure 1. In most current models 
block A is considered to be a single mode approximation to tank behavior - that is, a 
single mass (representing the fluid) together with appropriate springs coupling it to the 
remaining structure. However, it is clear that during a resonance of the fluid in the 
lower tank, lateral motion of the lower tank walls will influence the upper tank base mo- 
tion and hence, the upper tank base pressure. This type of relation is also required to 
properly couple the upper tank to the payload. 

This study presents a closed form approximation for the head to base response for 
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a rigid base cylindrical tank partially filled with fluid. This analysis considers the dis- 
tributed effects and the fluid pressure influences on shell motion. The model is not in- 
tended for detailed stability analysis, but as a tool in the early stages of vehicle design 
or for use in a rough cut at stability analysis. 

An experiment was performed to support the analysis. A cylindrical plastic tank 
containing water was longitudinally excited. Resulting data were compared with the 
analytical model. 


TRANSFER FUNCTION ANALYSIS 

The following basic assumptions are made in the analysis of the shell-fluid system 
to obtain the desired response function: 

(1) Shell and fluid are treated as one-dimensional in the longitudinal direction, and 
shell is distributed but has radial symmetry; that is, breathing and torsional modes 
are not considered. 

(2) Shell bending and lateral dynamic effects are neglected; that is, laterally, the 
shell is a membrane. 

(3) Fluid pressure response is assumed to depend only on tank base motion and is 
independent of wall motion except through acoustic velocity. 

(4) Ullage gas pressure is assumed to be constant. 

(5) Tank base and head are assumed rigid. 

The coordinate system for the shell-fluid system is shown in figure 2. (Symbols are 
defined in appendix A. ) The longitudinal displacement u of the shell is given by 

9^u £ 9u _ 9^u _ v_ dw I 

9x 2 * 31 ’ D St 2 " R 5x 
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(See, for example, refs. 10 and 11 with damping terms added.) The right side of this 
equation is a distributed forcing function which arises from the lateral motion of the 
wall. Also the Kj(9u/9t) term is a distributed mechanical damping. The term 
(y/R)(9w/9x) arises from the fluid pressures normal to the tank wall. These fluid 
pressures are determined using the equation (ref. 6, eq. (12)) 


i!p = c 2 + K* iE 


9x 


dr 


9t 


( 2 ) 


where 


s 2-flf B + 25 

g V Et j 


and 


A \ E t 


The Laplace transform of equation (1) with respect to time t with all initial con- 
ditions zero is 




_ v 9W 


(3) 


in terms of the Laplace variable s. The solution of equation (2) for fluid pressure to 
base velocity for the boundary condition (p(l , t) = Constant) is (ref. 7, eq. (38)) 


P(x, s) _ p s sinh Vc 2 s 2 + K*s (Z - x)] 

V(o, s) g |/ c 2 g 2 ~ j£* g cosh^c^s^ + K*s ij 


Now for the constant ullage pressure assumption, perturbations in fluid pressure do not 
influence longitudinal stress. Also when considering the shell as a membrane and 
neglecting lateral dynamics, the radial deflection is given by 
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O < X < l 


W(x, 


s) = _R? p( X) s) = _ / H 2 p \ sV(0, s) sin h [Vc 2 s 2 + K*s (l - x)] 
Et \g E T/ y c 2 s 2 + K*s cosh (Vc 2 s 2 + K*s 


( 5 ) 


Hence, 


v dW(x, S ) _ J 

R 3x 


gEr 


[cosh 

[V 

3 

! s 2 + K*s (l - x)_ 

^ cosh 

(v 

'c 2 s 2 + K*s l) 



o < x < Z 


( 6 ) 


0 Z < x < L 

This can also be expressed as 


v 3W(x, s)_ i/Rp sV(o, s)cosh 

k 

c 2 s 2 + K*s 

(l -x) 

R 3x ET S cosh ( 


2 s 2 + K*s l) 

1 


(7) 


where ^(x - l) is a unit step function. 

Since s is a parameter, equation (3) can be written as 


d 2 U 


m 


_ (k,s + — s 2 lU = H(x, s) 

2 V 1 D J 


dx 


( 8 ) 


Applying the Laplace transform with respect to x, St (where ft is the respective La- 
place variable) yields 

/a 2 <*r(p, s) - j&U(0 + , s) - ^ (0 + , s) - ^s + — - s 2 ^J<V(fa., s) s) (9) 

and letting dU(0 + , s)/dx = q(s) gives equation (9) as 


<*r(p, s) = 


XXP-, s) + /aU(0, s) + q(s) 

)» 2 -( k i s + ^= 2 ) 


( 10 ) 


The inverse Laplace transform with respect to f- of equation (10) is (details of calcu- 
lations shown in appendix B) 
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f y (coshCagX + a^) cosh|~a 1 (^ - x)J cosh [~o' 1 (x - l) coshCo^x - a^l) 


2 «2 “l + “2 


i(a 2 x - aqtfl 

“1- “2 J 


*1 “2 “1 " “2 

+ U(o, s)cosh(apx) + sinh( n„x) for 0 £ x < l 


-1- pxiui^u^^ 
a 2 


U(x, s) 


cosh(o? 2 x + 0 fjZ) coshj^o^x - Z)J coshja^x - Z)J cosh(or 2 x - a^l) 


(ID 


Y i 

pcosh(o? 2 x + afjZ) 

2a 2 

[_ “1 + “2 


*1 + a 2 


*1 - 2 


' a 2 


U(o, s)cosh(a 9 x) + sinh(d'px) for l < x < Z 


where 


<*1 = V 


2 2 T ,* 
cs + K s 


“ 2 =^ s 2 +K 1 ! 


_ yRp sV(o, s) 

ETg coshCajO 


At the tank top (x = L > £ ) 

fcosh(a-2L + a^l) 20^ coshfo^*- 1 - £)j cosMa^L - a^Z) 


U(L, s) = 


_ r 


2a 2 ^ a l + a 2 


2 2 
a l ~ a 2 


o-l - q-2 


+ U(o, s)cosh(o- 9 L) sinh(o- 9 L) for x = L (12) 


O' r 


The function q(s) required for these equations is determined by consideration of the 
head boundary conditions. 
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Head Boundary Condition 


The head and mechanical components above it are considered in terms of a me- 
chanical impedance Z m ; hence, 


z = F = — F — 

m V sU(L,s) 


(13) 


Then from the boundary condition at the top of the tank (ref. 7, for example) 


F = sZ m U(L, s) = - 27rRrE ( — - 


^ _ 2 \dx R / 


(14) 


x=L 


Since the radial deflection at x = L is zero and after rearranging 


dU 

dx 


= C 1 (s)U(L, s) 


(15) 


x=L 


where 


C x (s) = 


sZ md - * > 


2ttRtE 


(16) 


Now, evaluating dU/dx(L, s) and U(L, s) from equations (11) and (12), respectively, 
substituting into equation (15), and solving for q(S)/U(o, s) yield 


q(s) = 
U(o, s) C x 


! /3s 


2 Jsinf^OgL + a^l) 2a>2 sinhj^a^L - £)] 


a r 


sinh(«2L) - cost^a^L) 


«1 + a 2 


2 2 
<*l - a 2 




sinh(Q>2L - a-^l) 
a l - a 2 


Cips 

< 


|cosh(Q>2L + OjZ) 202 cosh [®2^ L " ^ 


2o 2 ^ «i + « 2 


2 2 
a 1 " a 2 


cosh(o2 L - l) 


a l~ <*2 


+ (»2 sinh(»2 L ) - Cj_ cosh^a-gL)] 


(17) 


I 



where for sinusoidal steady state V(o, s) = sU(o, s), and 


0 = 


yRp 


sV(0, s) s 2 d ^ s j Erg coshfo-jZ) 


(18) 


Transfer Function 


The desired transfer function can now be formed from equation (11) by dividing by 
U(0, s); that is. 


U(x, s) _ J 

U(o, s) 


0 g 2 fcost^o-gX + ct]l) 2«2 cosh^a-^Z - x)j cosh(a- 2 x - a-jZ )j 




2a- 0 I a^ + a 2 


2 J- 

a l ~ a 2 


, v sinh( a- 9 x) 
+ cosh(a- 9 x) + — SLJ 


a l ~ a 2 


for 0 < x < Z 


U(0, s) «2 


y (is) 


0 g 2 JcosM^gX + a-jZ) 2a?2 coshjj^Z - x)J cosh(a- 2 x - a-jZ)] 
2 £v 0 I o-j + ex £ „,2 _ „2 o-i - 


2 _2 
®1 - *2 


/ x sinh(o- 2 x) 
+ cosh(a 9 x) + f_ 


o-i- 0-2 


for Z < x < L 


U(0, s) CKg 




where q(s)/U(0, s) is defined by equation (17). For the tank head dynamics, the spe- 
cialized solution for x = L becomes 

U(L s) /Is** Jcosh(a- 2 L + a-jZ) 2a»2 coshja^L - Z)J coshfa^L - a-jZ) 

U(0, s) 2<* 2 arj + o- 2 


2 2 
o-l- a-2 


o-l - a-2 


/v sinh(a- 9 L) 
+ cosh(a- 9 L) +-2i^ _ — 


U(0, s) 0*2 


(20) 


This equation relates the tank head response to tank base input. Equation (19) and the 
special case equation (20) are the central result of the analysis. 
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Resonant Frequency Analysis 


It is of some interest to identify the resonant frequencies of the head motion rela- 
tive to the base. Considering the undamped case only (K* = 0, = 0) the resonant 

conditions can be found by the roots of the denominators of the various terms of equa- 
tion (20). 

Only two terms of equation (20) yield meaningful results. Setting the denominator of 
/3 (eq. (18)) equal to zero gives 


cosh ofjZ = cosh cs l = 0 


Replacing s with ico gives 


cos cc ol = 0 


which requires 


oj 


n 



or 





( 21 ) 


( 22 ) 


which are the modes associated with the fluid resonances (ref. 7). 

Considering the denominator of the term q(s)/U(0, s) sinl^argL)/*^] gives (where 
q(s)/U(o, s) is defined by eq. (17)) 


or 


C-, sinh(o' 9 L) 

cosh^gL) = 0 

a 2 


(23) 



(24) 
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Now, since Cj is a function of s and involves the mechanical impedance, the 
of a single rigid mass head will be considered (i. e. , Z m = Ms). Then, 

C . . Ms 2 (1 - , 2 ) . , c 2 
1 2 jtRtE * 


where 


c _ M(1 - /) 
2 2ttRtE 


These results substituted into equation (24) with s = ico give 


CgW sinh 


After some manipulation we obtain the result 


( K \ K ( K N 

inhl't/ wLi) =-»/ o>i coshl-t/ a; Li 

\v d / \V D ) 


tan Y = — 
Y 


where 



and 


n m sL _ Shell mass 

'“'Q = 

DCg Head mass 

Hence the mechanical resonances of the head are given by 


n 



case 


(25) 


(26) 


(27) 


(28) 


( 29 ) 
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where Y n are the solutions of equation (27). The mechanical resonances can be iden- 
tified by equation (29) and the solutions of equation (27) for Y n . Plots of these solu- 
tions are presented in figure 3. 

Several refinements can be made for this type of analysis. First, the effect of 
variation of acoustic velocity with distance can be accounted for. Acoustic velocity 
varies with distance because the end plates restrain the shell motion and thereby ef- 
fectively increase acoustic velocity in that area. When accounting for the variation in 
wave propagation time it is expected that the fluid resonant frequencies with variable 
acoustic velocity would be given by 


f n = 


n-i 

2 


f 


2 / c(x)dx 
'0 


The variation of c (local reciprocal acoustic velocity) with axial distance x can be 
expressed as 


c(x) = 



+ *m 

R 3p / 


(31) 


where the symbols are as defined earlier. 

The volume change due to shell bending can be calculated from classical shell 
theory (ref. 12, p. 36). For example, the solution for a cylindrical shell with the 
boundary conditions 


is 


x = 0 — w = 0 and ^ = 0 ^ 

dx 

x = L — w = 0 and — = 0 

dx 


(32) 


— dR _ 2R _ /2j// gin ginh + 2X cog ^ cos h <py] (33) 

R3p Et \R R / 


12 



where 




p 

2R /sin a cosh a - cos tv sinh 
Et V sin 2 or + sinh 2o> 


*) 


X = 


2R' 

Et V 


/cos a s inh a + sin a cosh a \ 
V sin 2cv + sinh 2 a / 


a 



y ( 3 « 


D 22 


y = x 


h 

2 


J 


Basically only the ends of the tank are influenced by this variation of c with length. 

A second refinement would be the inclusion of acoustic modes in the ullage gas. The 
frequencies of these modes are predicted by the equation 


f 


n, a 


na 

2(L - l) 


(35) 


where a is the acoustic velocity in the gas column and (L - l) is the gas column length. 
These gas column acoustic modes have a very narrow band width and when tuned with a 
fluid mode tend to depress it (ref. 13). It appeared that inclusion of this type mode in 
the analysis was inconsistent with the objectives of the present work. 


EXPERIMENTAL STUDY 

An experimental study was conducted to serve as a basis of comparison with the 
analysis. The test specimen was a cylindrical Plexiglas tank having a fixed rigid base 
and a removable rigid head. Plexiglas was chosen so that the fluid surface modes 
might be observed and for ease of fabrication. 

The salient characteristics of the tank are as follows: 
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Outside diameter, in. (m) 18 (0. 457) 

Wall thickness, in. (m) 1/8 (0.003175) 

Inside length, in. (m) 38. 75 (0. 984) 

Overall height, in. (m) 44.0(1.115) 

Weight of base and tube, lbm (kg) 86. 31 (39. 1) 

Weight of head, hardware, and instrumentation, lbm (kg) 28. 0 (12. 7) 

Density, lbm/ft^ (kg/m^) 74.3 (11 850) 

Material Plexiglas 


5 

The tank is shown mounted on a 28 000-pound force (1. 25x10 Newton) electrodynamic 
shaker in figure 4. Lateral support was not required with this shaker, nor was a sus- 
pension system. 

The instrumentation for the tests consisted of 11 piezoelectric accelerometers lo- 
cated to detect longitudinal motion (and 4 accelerometers for lateral motion of head and 
base). Also, two transducers measured tank base pressure at the center and one-half 
radius positions. 



Figure 4. - Plexiglas tank on shaker. 
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The normal testing procedure was to drive the system sinusoidally maintaining a 
1-g base acceleration amplitude (longitudinal) while sweeping from 20 to 400 hertz at a 
sweep rate of 1 octave per minute. The only variation from this procedure was for the 
mode shape studies where a selected constant resonant frequency was dwelled on for 
20 seconds. 

All test data were recorded on magnetic tape and analyzed using analog frequency 
analysis techniques. 


COMPARISON OF ANALYSIS AND EXPERIMENT 


The fundamental objective of the study was to predict the head to base transfer func- 
tion for a cylindrical tank partially filled with fluid. For this configuration the tank base 
pressure to base acceleration, and tank head acceleration to base acceleration transfer 
functions were analytically and experimentally determined for each 1/8 fullness point. 
Several parameters required in the analysis were either unknown or uncertain. The 
basic unknowns were the distributed damping numbers K* and and the reciprocal 
acoustic velocity of the fluid c. 


/-Local c(x) 



(a) Reciprocal acoustic velocity c as function 
of nondimensional fluid height III. 



(b) Damping number K* as function of non- 
dimensional fluid height III. 


Figure 5. - Parameters of analytical model as function of nondimensional 
fullness. 
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The difficulty in determining acoustic velocity is two-fold. The projected effects 
due to end boundary constraints were discussed previously. Secondly, uncertainty in 
selection of the proper modulus of elasticity from stress- strain diagrams was exper- 
ienced due to the plastic behavior of the Plexiglas. 

For these reasons the fluid reciprocal acoustic velocity c and the distributed 
damping number K* were chosen to approximate the fundamental fluid mode as iden- 
tified on the tank base pressure to tank base acceleration transfer functions. 

These parameters and a modulus of elasticity of 8. 6x10 pounds per square foot 
(4. 12x10^ N/m^) based on a mean c value (c = 0. 00181 sec/ft, 0. 00594 sec/m) were 
then used to predict the head to base transfer functions. For the distributed mechanical 

O O 

damping number a value of Kj = 0. 000015 second per square foot (0. 1615x10 sec/m ) 
approximated the results for all fullnesses. 

The c and K* values determined by the previous procedure have been plotted 
against relative fullness {l / L) in figure 5. The effective c of the system (fig. 5(a)) 
gradually decreases with increasing fluid height. Also shown on the plot is c(x) as pre- 
dicted by equations (31) to (34) and 


c 




c(x)dx 


In figures 6 to 10, respectively, typical results are presented for the 0, 1/4, 1/2, 
3/4, and 0. 96 fullness experimental and corresponding analytical transfer functions. 
The (a) and (b) figures are the head to base transfer function and the (c) and (d) figures 
are the base pressure to base acceleration transfer function. Mode identification for 
these plots will be aided by referring ahead to the summary plot figure 13. 

While the asymmetric modes were not predicted, the model did a reasonable job of 
predicting general response amplitude, indicating that the damping form is acceptable. 
The phase angle was apparently strongly influenced at the higher frequencies by the 
asymmetric modes. Apparent also in the head response data are the modes associated 
with the acoustics of the ullage gas. It will be noticed (i. e. , 1/4 fullness, fig. 7) that 
while these modes influence head motion the effects on tank base pressure are minimal. 

It is important to note that for some fullness conditions the contribution to the head 
response by the fluid mode is more important than that of the mechanical mode. This 
occurs for the quarter full case where the fluid mode (f = 152 Hz) has an amplitude 
ratio of 11 while the mechanical mode (f = 202 Hz) is about 5. 5. This also occurs for 
the 5/8 full case (data not presented) where the second fluid mode is dominant over the 
mechanical mode. This clearly indicates the importance of coupling the fluid behavior 
to the tank walls for predicting head motion. 
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Amplitude ratio, |U(x, s)/U(0, s)| 









Amplitude ratio, |p(0,s)/U(0, s)|, psi/g 









Heavy line Experiment 

Light line Analysis 


140 


100 



Frequency, Hz Frequency, Hz 


(c) Tank base pressure to tank base acceleration response; id) Tank base pressure to tank base acceleration response; 

amplitude ratio. phase angle. 

Figure IQ. - Concluded. 




It is also of interest to note for the 3/4 full case the coalescence of the second 
fluid mode and the head mechanical mode for the experiment while two separate res- 
onances showed for the model. The model resonance would coalesce at a slightly less 
full condition or with a small shift in acoustic velocity. 

A comparison of the fluid acoustic velocity assumptions is made in figure 11. Here, 
the fundamental fluid mode frequency is plotted against nondimensional fluid height for 
the constant and variable reciprocal acoustic velocity cases and the experiment. Either 
approximation is quite good for 1 / L > 0. 35. At this fluid height the wave length, 



Figure 1L - Comparison of experimental and analytical fundamental 
fluid mode frequencies. 


based on c = 0. 00181 second per foot (0. 00594 sec/m), is equal to the tank circumfer- 
ence; under such conditions the acoustic scattering becomes significant (ref. 4, p. 350), 
distorting and interfering with plane wave propagation. 

For an acoustic velocity of 1135 feet per second (346 m/sec) in the air above the 
liquid the frequencies of the gas column acoustic modes as predicted by equation (35) 
are plotted in figure 12 along with the experimentally observed modes. 

A summary of the experimentally observed resonant frequencies as a function of 
fullness is presented in figure 13. Identified on this plot are the fluid modes, the me- 
chanical head modes, and the ullage gas modes. 

Typical experimental modes shape plots generated by driving at the resonant fre- 
quencies are presented in figure 14 for the 3/4 full case. Here, shell longitudinal 


26 






Phase angle, Amplitude ratio, |U(x,s)/U(0,s)| Phase angle, Amplitude ratio, 

, s) to 11(0, s), deg U(x, s) to U(0,s), deg |l)(x, s)/U(0,s)| 



(a) Fluid fundamental mode. 



Nondimensional axial distance, x/L 


(b) Mechanical mode. 

Figure 14. - Shell longitudinal mode shape for 3/4 full tank. 




Figure 15. - Analytical pressure to base acceleration mode shape for 3/4 
full tank. 


acceleration at a given x-position to base acceleration, amplitude and phase is plotted 
against tank length. Also presented are the corresponding analytical results. The 
agreement between the two is good. The pressure to base acceleration mode shape was 
not measured in the experiment; however, the analytical values for these are presented 
in figure 15. 


CONCLUSIONS 

1. A model has been developed which predicts the head to base transfer function for 
a cylindrical tank partially filled with fluid. The model yields closed form results and 
accounts for the influence of fluid pressure effects and distributed mechanical and fluid 
damping. 
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2. The model did an adequate job of matching the transfer functions and mode 
shapes consistent with the objectives of the study. 

3. There remains some question as to proper choice of reciprocal acoustic velocity 
c and its variation with the x coordinate and fluid height. However, constant c values 
should be satisfactory. Also the distributed damping numbers and K* must be 
determined. 

4. The results clearly demonstrate the importance of considering fluid resonance 
effects on head to base response. Indeed in at least two cases the fluid mode effects 
dominated the mechanical modes for this transfer function. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 16, 1970, 

124-09. 


30 



APPENDIX A 


SYMBOLS 

A area, ft 2 (m 2 ) 

a acoustic velocity, ft/sec (m/sec) 

B fluid compressibility, ft 2 /lb (m 2 /N) 

c reciprocal acoustic velocity, sec/ft (sec/m) 

Cj function of s, see eq. (25) 

Cg constant, see eq. (26) 

Cg constant, see eq. (28) 

D rigidity of the shell, Er/(1 - y 2 ), lb/ft (N/m) 

E modulus of elasticity, lb/ft 2 (N/m 2 ) 

F force, lb (N) 

f frequency, Hz 

9 9 

g gravitational constant, ft/sec (m/sec ) 

H function defined by eq. (7) 

K proportionality constant per unit length for fluid damping, lb-sec/ft 2 (N-sec/m 2 ) 
K* distributed fluid damping, (K/A)[b + (2R/E r)], sec/ft 2 
Kj distributed mechanical damping, sec/ft (sec/m ) 

L tank height, ft (m) 

<£ Laplace transform operator with respect to t 

x Laplace transform operator with respect to x 

l height of fluid in tank, ft (m) 

M head mass, slugs (kg) 

m g mass per unit area of shell, p w t/s, lb-sec 2 /ft 3 (N-sec 2 /m 3 ) 
n mode number, integer 

p pressure, lb/ft 2 (N/m 2 ) 

p, Laplace variable with respect to distance x, 1/ft (1/m) 
q dU(0 + , s)/dx 

R radius, ft (m) 
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1 

$ 

unit step function 

s 

Laplace variable with respect to time, 1/sec 

t 

time, sec 

u 

displacement, ft (m) 

V 

velocity, ft/sec (m/sec) 

w 

radial displacement, ft (m) 

X 

distance coordinate, ft (m) 

Y 

-^m g /D Leo, rad 

y 

distance coordinate from center of tank, ft (m) 

m 

mechanical impedance, lb-sec/ft (N-sec/m) 

Z l> Z 2 

dummy variables for hyperbolic identify 

cv 

<pL/2 

®1 

Vc 2 s 2 + K*s, 1/ft (1/m) 

cv 2 

y(m g /D)s 2 + K x s, 1/ft (1/m) 

/3 

y/sv(0, s), see eq. (18), sec 2 /ft 2 (sec 2 /m 2 ) 

r 

(pRv/ETg) [sV(0, sJ/cosMcv^ )] , 1/ft (1/m) 

A 

variable defined by eq. (34) 

V 

Poisson ratio 

h 

dummy variable 

P 

density, lb/ft 3 (kg/m 3 ) 

T 

wall thickness, ft (m) 

<P 

- A l/«(l/m) 


V R 2 T 2 

* 

variable defined by eq. (34) 

CO 

angular frequency, rad/sec 

Subscripts: 

a 

acoustic mode 

n 

mode number, integer 

ef 

effective 
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x with respect to x 
w wall 

For variables U, P, H, V, and W, lower case represents the untransformed var- 
iable. Upper case represents variable transformed with respect to t, and script rep- 
resents variable transformed with respect to both t and x. 
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APPENDIX B 


INVERSE TRANSFORM CALCULATION 

To find the inverse Laplace transform with respect to x of equation (10), 


♦ S ) s) +q(s) 


p2 '( KlS+ TT s2 ) 


s) will be considered in three parts; namely, 


s) = 


s) 


m 


f* - 1*1 


s „2 


- K 1 s + — - s 


D 


(Bl) 


« 2 (p,s)= 


2 /„ m s 2) 

1 - k i s + F s 


(B2) 


and 


* 3 (p, s) = 


q(s) 


^ " ( K l s + ^ s2 J 


(B3) 


Equation (7) can be rewritten as 


H(x, s) = H 1 (x, s) |^1 -^(x - O] 


(B4) 


where 


HjCx, s) = 


vRp\ sV(0, s) 


ETg ' cosh (^c 2 s 2 + K*s z) 


cosh 


[Vc 2 s 2 + K*s (l - x)] (B5) 
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The inverse of ^(p-, s) then can be expressed as 


U x (x, s) = .sr 


Ji^H^x, s)[l - J(x - i)] 


2 ! v s 2 

? -( K l s+ — 8 


Making the substitutions 


o-j = Vc 2 s 2 + K*s 


«2 =\T ^ + K ! S 


y - vRp sV(o, s) 
ETg cosh( 0 ;^) 


and applying the convolution theorem to equation (B6) give 


^(x, s) = f y cosh [“l<* - flJt 1 - - oj — sinh|^Qf 2 (x - |)Jd| 


Uj(x, s) = 


2 — coshJaj(Z - ^smb^o-gCx - |)Jd£ 


for 0 < x < l 


2— coshja^ - |)Jsinh|o'2(x - £)Jd| 


for x > l 


Using the hyperbolic identity 
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sinh Zj cosh Zg = - sinl^Zj + Zg) + - sinh(Z 1 - Z 2 ) 
2 2 


and performing the required integrations give 


y r cosh|a^(2 - x)J cosh(a 2 x + a^l) cosh ^(x - f )J coshfo-gX - o' 1 Z)'j 


2q, 2 “1 + “2 


for 0<s<! 


*1 + “2 


*1 ' “2 


U 1 (x, s) =< 


(B9) 


2a<, 


^ coshja 2 (x - f)J cosh(a 2 x + a>j£) cosh£ar 2 (x - i)] coshCagX - a 1 x)'j 


for x a l 


*1 + “2 


*1 + “2 


*1 - “2 


J 


The inverse transform of equation (B2) with respect to x is readily found as 

U 2 (x, s) = U(0, s)cosh(o! 2 x) (BIO) 

To obtain the inverse transform of equation (B3) it is noted that the unknown q is only 
a function of s; hence, 


Uq(x, s) = q(s) — sinh(a? 2 x) (Bll) 

a 2 

Combining the results of equations (B9), (BIO), and (Bll) gives the inverse transform 
of equation (10) as 


r , 

"cosh(of 2 x + otyf') cosh 

- x) 

cosh 

^(x - l) 

coshC^x - 

2 “ 2 ! 

“1 + “2 “1 + a 2 

a 

1 ■ “2 

“1 - «2 J 


+ U(0, s)cosh(a 2 x) + SM sinh( a^x) 
“2 


for 0 < x < l 


u(x, s) = *\ 


r J 

fcoshta'gX + a^l) cosh a 2 (x - £)j 

coshj 

a 2 (x - 0 ] 

cosh(of 2 x - a^l) | 

2 “2 

[ “1 + “2 “l + “2 

“1 - a 2 

“1 - “2 J 


f (B12) 


+ U(0, s)cosh(a2 x ) + Sinh( ff 2 x) 

a 2 


for L < x < L J 
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where 


c 2 s 2 + K*s 


« 2=\^ s2 + K 1 S 


y= vBp sV(o, s) 
ETg COShCo!^) 


At the tank top (x = L > l) 


U(L, s) 


y jcosh(o! 2 L + QfjO 2 q !2 cosh ^(L " Oj cosMa^L - 


a l + a 2 


2 2 
a l ~ a 2 


*1 - a 2 


+ U(o, sJcoshCagL) + sinh( L) for x = L (B13) 
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